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superlattices,  (ii)  interaction  of  hydrogen  with  metals  and 
hydrogen  absorption  around  metal  ions,  (iii)  defects  and  defect 
complexes,  (iv)  electronic  structure  of  LiMgH3  class  of  compounds 
and  (v)  super-shell  structure  in  clusters.  In  the  following  we 
provide  a  brief  outline  of  our  accomplishments  in  each  of  these 
areas  as  well  as  our  professional  accomplishments. 


B.  SUMMARY  OF  THE  IMPORTANT  RESULTS. 

(i)  Magnetism  of  Low  Dimensional  Systems.  Transition 

Metal  Superlattices  and  Clusters. 

Considerable  effort  was  devoted  to  understanding  the  effect  of 
size,  dimensionality  and  temperature  on  the  magnetic  moment  of 
transition  metal  systems1,2,7,9,27.  These  include  Fe,  Co  and  Ni  forming 
chains,  planes,  slabs  of  varying  thickness,  clusters,  and  bulk 
containing  vacancies,  impurities  and  vacancy-impurity  complexes. 
Since  the  existing  theories  are  designed  to  study  periodic  systems, 
we  developed  a  new  method  by  combining  ab-initio  and  real-space 
tight  binding  methods1,2.  We  found  that  the  magnetic  moment  of  an 
atom  is  primarily  governed  by  the  local  coordination  and  that 
chains  are  more  magnetic  than  planes  which  are  more  magnetic  than 
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bulk.  Our  studies  on  chains  of  finite  number12  of  atoms  showed  that 
the  Ni  chains  containing  odd  number  of  atoms  had  higher  moments 
than  even  number  chains.  This  quantum  size  effect12  persists  in 
chains  having  up  to  20  atoms  at  which  point  the  moments  approach 
the  asymptotic  limit.  The  appearance  of  a  non  bonding  orbital  at 
the  Fermi  energy  was  found  to  account  for  the  larger  moment  in  odd 
atom  chains. 

Calculations  were  carried  out  on  magnetic  clusters  modelled  as 
fragments  of  the  bulk2.  It  was  found  that  for  ferromagnetic 
materials  like  Fe,  Co,  or  Ni,  the  clusters  have  higher  moment  than 
bulk  and  the  moment  on  an  atom  approaches  the  bulk  limit  when  an 
atom  is  surrounded  by  three  neighboring  shells.  We  also  found  that 
mono-layers  and  small  clusters  of  non  magnetic  elements  are 
magnetic  like  V  are  magnetic19. 

Stern-Gerlach  experiments  on  Fe„  and  Con  clusters  in  beams 
showed  that  these  clusters  have  much  lower  magnetic  moments  per 
atom  than  the  bulk.  This  observation  defied  the  theoretical 
predictions  of  higher  magnetic  moment  in  smaller  sizes.  We 
resolved  this  apparent  paradox18,22  by  showing  that  small  clusters 
behave  like  paramagnetic  atoms  and  that  the  observed  small  moment 
is  a  consequence  of  the  superparamagentic  relaxation  at  the  cluster 
temperatures.  Our  model  explained  the  experimental  data  on  the 
size  and  temperature  dependence  of  the  magnetic  moment  and  has 
since  been  used  by  various  experimental  groups  to  explain  their 
experiments.  We  also  studied  the  magnetic  behavior  of  small 
clusters  with  spins  coupled  antiferromagnetically26'28.  Here  we 
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showed  that  the  geometry  had  a  strong  effect  on  the  observed 
behavior.  While  cuboctahedric  clusters  had  well  defined  ground 
states,  the  icosahedric  clusters  have  frustrated  free  spins  due  to 
the  underlying  lattice  and  hence  will  behave  as  paramagnets. 

We  have  also  studied  the  magnetism  in  transition  metal 
super  lattices29.  Experimentally  the  Fe  layers  separated  by  Cr  space 
layers  have  been  found  to  couple  antiferromagnetically .  To 
understand  this  finding,  we  have  studied  the  magnetic  coupliny 
between  Fe  layers  separated  by  spacer  layers  consisting  of  3d 
transition  metal  elements  (Sc,  Ti,  V,  Cr,  Fe,  Co,  Ni  and  Cu)  by 
using  the  self-consistent  density  functional  theory  and  discrete 
variational  method.  The  coupling  between  the  Fe  layers  is 
ferromagnetic  in  all  cases  except  for  Cr  where  this  coupling 
alternates  from  ferro  magnetic  to  antiferromagnetic  depending  on 
whether  the  spacer  layers  are  odd  or  even.  Furthermore  the  spacer 
layers  involving  Sc,  Ti  and  V  are  antiferromagnetically  coupled  to 
Fe  while  Co  and  Ni  layers  are  coupled  ferromagnetically .  These 
findings  are  the  subject  of  experimental  investigations  in  various 
labs. 

(ii)  Interaction  of  Hydrogen  with  Metals  and  Hydrogen 

Absorption  Around  Metal  Ions 

In  most  metals  hydrogen  precipitates  at  low  temperatures  to 
form  an  ordered  hydride  phase.  A  departure  from  this  common 
occurrence  was  first  observed  in  resistivity  experiments  in  rare- 
earth  metals  (Sc,  Y,  etc.).  The  ability  of  these  materials  to 
retain  hydrogen  in  solid  solution  down  to  absolute  zero  was 
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completely  unexpected  and  gave  rise  to  a  flurry  of  experimental 
studies  involving  neutron-scattering,  nuclear  magnetic  resonance 
and  channelling.  The  above  anomaly  was  shown  to  be  due  to  the 
pairing  of  hydrogen,  although  the  structural  origin  of  the  pairs 
remained  a  debatable  issue. 

Using  clusters  of  Y  as  a  model  of  the  bulk,  we  calculated  for 
the  first  time  self-consistently  the  energetics314  of  a  pair  of 
hydrogen  atoms  located  at  adjoining  tetrahedral  (T)  and  octahedral 
(0)  sites,  T-T  sites  and  T  sites  mediated  by  a  metal  atom.  The 
last  configuration  was  found  to  be  energetically  most  stable  in 
agreement  with  all  available  experimental  data.  Subsequent 
calculations  have  confirmed  our  results.  The  binding  energy  of  the 
pair  thus  calculated  was  used  in  a  lattice  gas  model620  to  predict 
the  temperature  dependence  of  the  fraction  of  H  atoms  involved  in 
pairs.  This  prediction  has  now  been  verified  in  inelastic  neutron 
scattering  experiments. 

Hydrogen  absorption  in  metals  is  extremely  important  for 
industrial  applications.  We  have  carried  out  Quantum  chemical 
calculations25  including  correlations  to  study  the  amount  and  nature 
of  hydrogen  uptake  by  a  Li  atom  in  neutral  as  well  as  in  ionic 
form.  We  have  found  that  while  the  neutral  Li  atom  cannot  bind 
more  than  one  hydrogen  atom,  the  cation  can  bind  a  large  number  of 
hydrogen  atoms  in  molecular  form.  This  contrasting  behavior 
originates  from  the  fact  that  the  bonding  of  hydrogen  with  neutral 
Li  is  covalent  while  that  with  the  cation  is  primarily  caused  by 
polarization.  The  bonding  of  hydrogen  with  the  Li  ion  exhibits 
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dual  nature.  The  first  H2  molecule  binds  to  the  anion 
dissociatively  while  the  subsequent  H2  molecule  can  only  exist  in 
a  "physisorbed"  state.  Our  findings  have  just  been  confirmed  by 
Prof.  Bowers  at  UCLA  who  has  been  able  to  attach  up  to  7H2  molecules 
to  a  0o+  ion. 

(iii)  Defects  and  Defect  Complexes  in  Metals: 

Using  the  density  functional  theory  and  the  embedded  cluster 
model,  we  addressed  a  long  standing  problem  of  the  effect  of  carbon 
impurities  on  the  migration  and  trapping  of  vacancies  in  Fe.  We 
calculated11,  self-consistently  for  the  first  time,  the  energetics 
of  carbon-vacancy  complexes  in  Fe  as  a  function  of  distance 
separating  the  carbon  atoms  and  the  vacancy  along  different 
crystallographic  directions.  Carbon  was  found  to  prefer  an  off- 
center  site  from  the  vacancy  irrespective  of  its  direction  from  the 
vacancy  center.  The  equilibrium  site  of  carbon  was  one  where  it 
was  coordinated  to  three  Fe  atoms  in  its  nearest  neighbor  shell. 
Multiple  decoration  of  vacancies  by  carbon  was  also  found  to  be 
possible.  Our  results  were  consistent  with  positron  annihilation 
experiments.  We  argued  that  positron  annihilation  spectroscopy 
cannot  provide  the  geometry  of  the  defect  complex  and  that 
experiments  more  sensitive  to  the  geometry  of  the  carbon-vacancy 
complexes  are  needed  to  compare  with  our  theory.  We  showed  that 
the  bonding  of  carbon  to  iron  atoms  was  not  caused  by  charge 
transfer  between  the  atoms,  but  rather  had  a  magnetic  origin. 

We  were  the  first  to  resolve  an  apparent  conflict6  between  the 
perturbed  angular  correlation  (PAC)  experiment  and  the  kinetic 
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model  concerning  the  equilibrium  site  of  In  on  Cu  (100)  and  the 
clustering  of  In  atoms  at  elevated  temperatures.  A  comparison  of 
the  calculated  and  experimental  result  for  the  electric  field 
gradient  (EFG)  and  asymmetry  parameter  (t))  led  to  the  conclusion 
that  In  occupies  the  substitutional  terrace  site  on  Cu  (100) 
surface.  The  migration  and  clustering  of  In  atoms  at  elevated 
temperatures  obtained  from  the  same  experiment  was  explained  in  a 
kinetic  model  by  assuming  In  atoms  to  lie  on  a  surface  that  has 
over  90%  vacant  sites.  Clearly  if  In  occupied  the  substitutional 
terrace  site,  other  sites  on  the  same  surface  are  not  empty  but 
occupied  by  Cu-atoms.  Thus,  there  was  a  conflict  between  the 
assumed  equilibrium  In  site  and  the  nature  of  diffusion  of  In  at 
higher  temperatures.  This  conflict  was  resolved5  by  calculating  the 
energies  of  In  occupying  the  on-top  site,  bridge  site  and  the 
substitutional  terrace  site.  In  each  case  the  distance  of  In  from 
the  substrate  was  varied  to  minimize  the  energy.  We  found  that  the 
bridge  and  substitutional  sites  compete  for  In  trapping,  with  the 
bridge  site  being  slightly  favored  over  the  substitutional  site. 
As  the  sample  is  annealed,  some  of  the  In  atoms  in  bridge  sites 
migrate  to  the  substitutional  terrace  site  where  efg  is  observed. 
The  clustering  arises  as  the  In  atoms  at  the  bridge  sites  move 
freely  and  interact  with  the  In  atoms  at  the  terrace  site. 
Subsequent  to  our  calculation,  the  bridge  site  for  In  was 
identified  experimentally. 

Encouraged  by  the  success  of  the  cluster  model  in  depicting 
the  electronic  structure  and  binding  at  surfaces,  we  turned  our 
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attention  to  understanding  the  nature  of  bonding*16  at  the  interface 
of  metal  matrix  composites  (Al/SiC) .  There  has  been  considerable 
debate  concerning  whether  the  bonding  is  due  to  charge  transfer  or 
diffusion  of  atoms  across  the  interface.  Recently,  experiments 
involving  Auger  and  electron  loss  spectroscopies  have  been  used  to 
determine  the  mechanism  of  formation  of  the  Al/SiC  interface. 
Annealing  at  moderate  temperatures  was  found  to  cause  aggregation 
of  A1  at  C-rich  sites.  Furthermore,  a  shift  of  the  A1  2p  level  to 
higher  binding  energy  and  Si  2p  level  to  lower  binding  energy  was 
also  observed.  We  have  demonstrated16  through  SCF-LCAO-MO 
calculation  that  all  of  these  results  can  be  understood  by  studying 
energies,  electron  charge  transfer,  bond  lengths,  and  core-level 
shifts  of  Al-Al,  Al-C,  Al-Si  and  Si-C  dimers.  Our  results  also 
revealed  that  the  bonding  at  the  interface  can  be  explained  by  a 
charge  transfer  mechanism. 

Disordered  materials  constitute  and  extreme  form  of  defects. 
We  have  looked17  at  the  electronic  conduction  in  disordered  metals 
and  alloys.  We  developed  and  electronic  transport  scheme  which 
could  permit  dc  and  ac  conductivity  calculations,  could  be  applied 
to  the  localized  and  extended  regimes,  and  could  be  combined  with 
ab-initio  tight  binding  scheme  to  allow  ab-initio  transport 
calculations.  We  used  the  Kubo-  Greenwood  formula  for  conductivity 
and  the  theory  of  orthogonal  polynomial.  The  new  formulation, 
developed  within  the  tight-binding  model,  provides  a  systematic 
expansion  beyond  the  random  phase  approximation  and  permits  studies 
of  ac,  dc  conductivities,  magnetoresistance  and  the  correlation 
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functions  (velocity-velocity  correlation  function) .  Since  the 
formulation  is  carried  out  in  real  space,  it  can  be  applied  to 
disordered  systems  with  ease.  We  have  used  our  approach  to  study 
the  canonical  site-diagonal  disordered  cubic  lattice,  and  the 
results  are  in  perfect  agreement  with  other  approaches.  We  are  now 
planning  to  use  it  to  perform  studies  on  clusters  and 
quasicrystals . 

(iv)  Electronic  Structure  of  LiMqH^  Class  of  Compounds 

The  possibility  of  LiMgH3  class  of  materials  showing 
superconductivity  at  elevated  temperatures  has  attracted 
considerable  interest.  We  have  calculated23  the  lattice  spacing  and 
the  electronic  structure  of  LiMgH3,  NaMgH3  and  LiCaH3  by  modelling 
these  hydrides  as  clusters  with  modified  perovskite  structure  and 
using  the  self-consistent  field  linear  combination  of  atomic 
orbitals  —  molecular  orbitals  theory.  This  yields  lattice 
constants  of  3.47  A,  3.53  A  and  4.07  A  for  LiMgH3  and  LiCaH3 
respectively.  The  electronic  structure  of  LiMgH3  is  characteristic 
of  an  insulating  material.  The  gap  at  the  Fermi  energy  in  the 
density  of  states  of  LiMgH3  is,  however,  significantly  less  than 
that  of  LiBeH3  and  is  insensitive  to  moderate  changes  in  the 
interatomic  spacing.  On  the  other  hand,  LiCaH3  and  NaMgH3  are  both 
found  to  exhibit  metallic  character  with  the  Fermi  energy  lying  in 
regions  of  high  electron  density  of  states. 

(v)  Super-Shell  Structure  in  Clusters 

We  examined  the  possibility  that  clusters  may  exhibit  super¬ 
shell  structures21  with  orbital  angular  momentum,  i> 6,  that  has  long 
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been  predicted  for  nuclei,  tve  have  estimated  an  upper  limit  for 
the  number  of  atoms  by  comparing  the  energy  difference  between  the 
highest  occupied  and  the  lowest  unoccupied  state  with  the  crystal 
field  splitting  of  the  cluster.  The  former  was  obtained  by  solving 
the  Schrodinger  equation  for  a  spherical  potential  well  with  hard 
walls,  while  the  latter  is  obtained  from  the  band  structure  of  the 
solid.  The  results  indicate  that  shell  structure  may  persist  in 
clusters  containing  as  many  as  a  million  atoms.  Experimentally  the 
super-shell  structure  has  already  been  seen  in  clusters  of  500 
atoms . 
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Stockholm,  Sweden  1990 

International  Symposium  on  "Metal-Hydrogen 

Systems",  Banff,  Canada  1990 

Chinese  Materials  Research  Society  Meeting, 

Beijing  (declined)  1990 
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